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Abstract. We consider stochastic differential equations dY = V (Y) dX driven 
by a multidimensional Gaussian process X in the rough path sense. Using 
Malliavin Calculus we show that Yt admits a density for t £ (0, T] provided 
(i) the vector fields V = (Vi, Vd) satisfy Hormander's condition and (ii) the 
Gaussian driving signal X satisfies certain conditions. Examples of driving sig- 
nals include fractional Brownian motion with Hurst parameter H > 1/4, the 
Brownian Bridge returning to zero after time T and the Ornstein-Uhlenbeck 
process. 



1. Introduction 

In the theory of stochastic processes, Hormander's theorem on hypoellipticity 
of degenerate partial differential equations has always been an important tool to 
decide whether or not a diffusion process with given generator admits a density. 
This dependence on PDE theory was removed when P. Malliavin devised a purely 
probabilistic approach to Hormander's theorem which is perfectly adapted to prove 
existence and smoothness of densities for diffusions given as strong solution to an 
ltd stochastic differential equation driven by Brownian motion. 

The key ingredients of Malliavin's machinery, better known as Malliavin Cal- 
culus or stochastic calculus of variations can be formulated in the setting of an 
abstract Wiener space (W, Tt, fi). This concept is standard (e.g. [3D] or any modern 
book on stochastic analysis) as is the notion of weakly non-degenerate K e -valued 
functional ip which has the desirable property that the image measure tp^fi, is ab- 
solutely continuous with respect to Lebesgue measure on R e . (Functionals which 
are non-degenerate have a smooth density.) Precise definitions are given later on 
in the text. 

Given these abstract tools, we turn to the standard Wiener space C ([0, T] ,R d ) 
equipped with Wiener measure i.e. the standard model for Brownian motion B = 
B (uj). From Ito's theory, we know how to solve the stochastic differential equation 

d 

dY = ^ V (Y) o dB l = V (Y) o dB, Y (0) = y e K e . 

i=i 

The Ito-map B t— > Y is notorious for its lack of strong regularity properties. On 
the positive side, it is smooth in a weak Sobolev type sense ("smooth in Malliavin's 
sense") and under Hormander's condition at ijq € IR e 

(1.1) (H): Ue[V u ...,V d } Vo =T yo R e =^ e 

one can show (e.g. (30J [36l El H2] ) that the solution map B i— ► Yt is non-degenerate 
for all t £ (0, T]. This line of reasoning provides a direct probabilistic approach to 
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the study of transition densities of Y and has found applications from stochastic 
fluid dynamics to interest rate theory, e.g. [T21ET]. The same range of application^ 
nowadays demand stochastic models of type 

(1.2) dY = V{Y)dX 

where X is a Gaussian process, such as fractional Brownian motion (short: fBm). 
Differential equations of this type have also been used as simple examples for the 
study of ergodicity of non-Markovian systemems, [24] , 

In a previous paper [4] we linked rough paths and Malliavin calculus by giving 
a (simple) proof that existence of a density for solutions of (II. 2p holds true under 
ellipticity i.e. 

(E) : Span [V 1 , V d ] yo = T yo R a = M e . 

and generic non-degeneracy conditions on X, the differential equation (II. 2|) being 
understood in the rough path sense [25l [29] , a unified framework which covers at 
once Young and Stratonovich solutions (and goes well beyond). The aim of this 
paper is to prove the existence of densities under Hormander's condition (H) in the 
following form: 

Theorem 1. Let (Xj, . . . , X t d ) = (X t : t <G [0, T]) be a continuous, centered Gauss- 
ian process with independent components X 1 , . . . , X d . Assume X satisfies the con- 
ditions listed in section [7} (In particular, X is assumed to lift to a geometric 
rough path so that i f _?.!?)) makes sense as random rough differential equations.) Let 
V = (Vi, ...,Vd) a collection of smooth bounded vector fields on R e with bounded 
derivatives which satisfies Hormander's condition (H) at Dq. Then the random 
RDE solution Y t = Y t (u) S R e to il.ty) started at Yq — yo admits a density with 
respect to Lebesgue measure on R e for all times t € (0,T]. 

One should note that X, the Gaussian driving signal of (|1.2[) , is fully described 
by the covariance function of each component and, under the further assumption of 
IID components, by the covariance of a single component, i.e. R (s, t) = E (X]Xl^ . 
In principle all conditions on X are checkable from the covariance, in practice 
it is convenient to have conditions available which involve the reproducing kernel 
Hilbert ~ or Cameron- Martin space associated to X as well as certain sample 
path properties. Leaving these technical details to section 0] we emphasize that 
our conditions are readily checked in many cases and now give a list of examples 
to which our theorem applies. It may be helpful to note that whenever X is a 
semi-martingale on [0, T] then (|1.2jl can be understood as Stratonovich stochastic 
differential equation, i.e. 

d 

dY = ^Vi{Y)odX l . 

i=l 

In such cases, rough path theory appears as intermediate tool that is neither needed 
to understand the assumptions nor the conclusions of theorem [TJ There may be 
cases when X can be written in terms of Brownian motion so that ultimately the 
techniques of [5, 38 j are applicable. But in general theorem Q] covers new grounds. 



^For instance, stochastic differential equations driven by fBM have applications applications 
to vortex filaments; applications to finance include (geometric) fractional Brownian motion as 
paradigm of a non-semimartingale which admits no arbitrarge under transaction costs. The reader 
is refered to the books 32] Sec. 5.3 and 5.4], [9| Sec 8.1.] and the references therein. 



DENSITIES FOR RDES UNDER HORMANDER'S CONDITION 



Example 1 (Brownian motion). When R (s, t) = min (s, t) the driving signal X — 
X (uj) is a d- dimensional standard Brownian motion on [0, T]. As is well understood 
[29) one needs to add Levy's area process to obtain a geometric rough path (known in 
this context as Brownian rough path or Enhanced Brownian motion). A solution to 
1 1.2) in the rough path sense then precisely solves the stochastic differential equation 
in Stratonovich form 

d 

dY = ^2v i (Y)odB i . 

Subject to Hormander's condition (H), Theorem\]\ then shows that Y t = Y t (lo) has 
a density for t > which is of course well-known. 

Example 2 (Fractional Brownian motion). When 2R (s, t) = t 2H + s 2H — \t — s\ 2H 
the corresponding process is known as fractional Brownian motion B H = B H (uj) 
with Hurst parameter H 6 (0, 1). Popularized by [31] it generalizes Brownian mo- 
tion which corresponds to H = 1/2. For fixed H and when no confusion is possible 
we shall write B = {B 1 ,...,B d ) for d-dimensional fBm. When H > 1/2, Kol- 
mogorov's criterion shows that B has nice sample paths (more precisely, Holder 
continuous sample paths of exponent greater 1/2) which has the great advantage 
that \1.2) can be understood as integral equation for fixed u) based on Young inte- 
grals (i.e. limits of Riemann-Stieltjes sums). In this setting of nice sample paths, 
existence of a density was established in [33] assuming ellipticity. Using deter- 
ministic estimates for the Jacobian of the flow this density was then shown to be 
smooth [34] : building on the same estimates the Hormander case was obtained in 
PQ. For H < 1/2 the situation appears to be fundamentally different: first, in view 
of Brownian (and worse) sample path regularity one needs ltd or rough path ideas 
to make sense of hi. 13) for H < 1/2. Secondly, the proof of pQ does not extend 
to the rough path setting and relies somewhat delicately on specific properties of 
fractional Brownian motion. In any case, theorem^ shows that Y t — Y t (lu), solu- 
tion to the RDE driven by multidimensional fBM with Hurst parameter H > 1/4 
has a density for all positive times provided the vector fields satisfy Hormander's 
conditio"^. The novelity is of course the degenerate regime H < 1/2 with sample 
path regularity worse than Brownian motion. 

Example 3 (Brownian Bridge). Let B be a d-dimensional standard Brownian mo- 
tion. Define the Brownian bridge returning to zero at time T by 

Xj : =B t -^B T forte [0,T]. 

Equivalently, one can define X T via the covariance 

R T {s,t) = min (s, t) (1 — max (s, t) /T) . 



2 The estimates of |34| can be generalized [20] to (sharp) deterministic estimates on the Jacobian 
of RDE solutions giving L p -estimatcs on the flow of RDEs driven by fBM if and only if H > 
1/2. In particular, one sees that L p -estimates on the flow of Stratonovich SDEs (H = 1/2) are 
fundamentally probabilistic i.e. rely on cancellations in stochastic integration. At present, the 
question of how to obtain good integrability when H < 1/2 is open although one suspects that 
Gaussian isoperimetry will ultimately play a role. 

■^As is well understood |29| , for H < 1/4 fractional Brownian increments decorrelate too slowly 
for stochastic area to exist and so there is no meaningful lift of fBM with H < 1/4 to a geometric 
rough path. 
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Clearly \t—T = and trivially (take dY = dX) the conclusion of Theorem^ 
cannot hold; this behaviour is indeed ruled out by condition [3| in section [^J On 
the other hand, we may consider X T+£ restricted to [0, T] and in this case the 
conditions in section^ are readily verified. It is worth remarking that Z := X T+E 
stopped at time T is also a semi-martingale; for instance, by writing (Xj +e : t < T) 
as strong solution to an ltd differential equation with (well-behaved) drift (as long 
as t < T). The conclusion of theorem]^ can then be stated by saying that the 
unique Stratonovich solution to dY = Y)Vj (Y) o dZ l admits a density for all times 
t £ (0,T] provided the vector fields satisfy Hormander's condition (H). 

Example 4 (Ornstcin-Uhlcnbcck). Let B be a standard d-dimensional Brownian 
motion and define the centered Gaussian process X by Wiener-Ito integration, 

X\ = [ e- (t - r) dB l r with i = 1, d. 
Jo 

X satisfies the ltd differential equations, dX t — —X t dt + dB t and is also a semi- 
martingale. The conditions of section^ are readily checked (in essence one uses 
X t ~ B t at t —t 0+ and the absence of Brownian Bridge type degeneracy). The 
conclusion of theorem [7] can then be stated by saying that the unique Stratonovich 
solution to dY = ^ Vi (Y) o dX 1 admits a density for all positive times provided 
the vector fields satisfy Hormander's condition (H). 

Further examples (for instance, "fractional" version of the Brownian Bridge and 
Ornstein-Uhlenbeck process) are readily constructed. Generalizing examples [21 [4] 
one could consider Volterra processes [8], i.e. Gaussian process with representation 
X = J K (■, s) dB s , and derive sufficient conditions on the kernel K which imply 
those of section 2J Existence of a rough path lift of X aside, one would need non- 
degeneracy of K and certain scaling properites as t — > 0+ but we shall not pursue 
this here. (In any case, there are non- Volterra examples, such as the Brownian 
bridge returning to zero at (T + e), to which theorem [T] applies.) 

The proof of theorem [T] is based on the fact [1] that RDE solutions driven by 
Gaussian signals are " 7i-differentiable" i.e. differentiable in Cameron-Martin direc- 
tions. Existence of a density is then reduced to showing that the Malliavin covari- 
ance matrix is weakly non-degenerate. The standard proof of this (e.g. [30l 12) or 
[321 Sec 2.3.2]) is based on Blumenthal's 0-1 law and the Doob-Meyer decomposi- 
tion for semi-martingales. The main difficulty to overcome in the general Gaussian 
context of this paper is that the Doob-Meyer decomposition is not available and 
we manage to bypass its use by suitable small time developments for RDEs, ob- 
tained in [19] ; in conjunction with (Stroock-Varadhan type) support description for 
certain Gaussian rough paths (as conjectued by Ledoux et al. [55] and carried out 

independently in [HI 02], see a ^ so H anc ^ 

The crucial induction step - which explains the appearance of higher brackets - 
requires us to assume a "non-standard" Hormander condition which involves only 
iterated Lie-brackets contracted against certain tensors arising from free nilpotent 
Lie groups. Equivalence to the usual Hormander condition (H) is then established 
separately. 
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2. Preliminaries on ODE and RDEs 

2.1. Controlled ordinary differential equations. Consider the controlled or- 
dinary differential equations, driven by a smooth i? d -valued signal / = / (t) along 
sufficiently smooth and bounded vector fields V = (Vi, Vd), 

d 

(2.1) dy = V(y)df = J2 V i (V) f (*) dt > V (*o) = Vo £ K e . 

»=1 

We call U(^ t (yo) = yt the associated flow. Let J denote the Jacobian of U. It 
satisfies the ODE obtain by formal differentiation w.r.t. yo. More specifically, 

a ^ \ le U ^ to ^° + ea ) 
is a linear map from K e — > R e and we let j/<_ t (yo) denote the corresponding e x e 



£ = 



to 

matrix. It is immediate to see that 



pU (yo) 



dt 

where • denotes matrix multiplication and 



JLt (vo) 



±Mf(y,t) = ^v>(y)p. 

i—l 

Note that j/ . t — jl ■ jf . We can also consider Gateaux derivates in the 
driving signal and define 



DuU f - I —U f+eh 



e=0 



One sees that DfJj(^_ Q satisfies a linear ODE and the variation of constants formula 
leads to 

DhUU (y ) = fj2 J Ls (Vi (uLo)) dK- 
Finally, given a smooth vector field W a straight-forward computation gives 
(2-2) dJ&_ t (W (U^o) ) = j£_ t ( [Vi , W] (c//_ ) ) C • 

2.2. Rough differential equations. Following [551 I2H1 HZ] a geometric p-rough 
path x over R d is a continuous path on [0, T] with values in G^ , the step- 
[p] nilpotent group over W d , and of finite p- variation relative to the [17] Carnot- 
Caratheodory metric d on G^ (R d ), i.e. 

sup sup y^^d(x ti ,K U+1 ) p < oo. 

neN0<ii<...<t 7l <T~^ 

As in [6l [25] we view G^ (ffi. d ) as embedded in its enveloping tensor algebra i.e. 
g Ip] ( R d) c T W ( R d) := e^o,...^] (R d )® fc • 
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One can then think of x as a path x : [0, T] — > M. d enhanced with its iterated 
integrals although the later need not make classical senstQ. The canonical projec- 
tion to (R d )® fc is denoted iik (x) or x fe . Lyons' theory of rough paths then gives 
deterministic meaning to the rough differential equation (short: RDE) 

(2.3) dy = V(y)dx. 

(One can think of RDE solutions as limit points of corresponding ODEs of form 
(12. ip in which the smooth driving signals plus their iterated integrals up to order 
[p] converge to x in suitable p- variation distance.) The motivating example, e.g. 
[25} 129]. is that almost every continuous joint realization of Brownian motion and 
Levy's area process (equivalently: iterated Stratonovich integrals) gives rise to a 
geometric p- rough path for p > 2, known as Brownian rough path or Enhanced 
Brownian motion (cf. example [lj which provides in particular a robust path-by- 
path view of Stratonovich SDEs. 

Back to the deterministic RDE (|2.3p and assuming smoothness of the vector 
fields V = (Vi,..., Vd), the solution induces a flow yo i— > Uf^ tf) (yo). Following 
[2Bl[17]i the Jacobian J£_ to of the flow exists and satisfies a linear RDE, as does 
the directional derivative 

Dull* - / —U T ' hX \ 
UhUt -°-\de Ut -° J £=0 

for a smooth path h. If x arises from a smooth path x together with its iterated 
integrals the translated rough path T^x (cf. (2S1 [29]) is nothing but x + h together 
with its iterated integrals. In the general case, we assume h € (7<?- var with l/p + 
1/q > 1, the translation X/jX can be written in terms of x and cross-integrals 
between 7Ti (x 0j .) —: x and the perturbation h. (These integrals are well-defined 
Young- integrals . ) 

Proposition 1. Let~K be a geometric p-rough paths overW 1 and h e C q - Var ([0,T] ,R d ) 
such that 1/p + 1/q > 1 . Then 

D h U* (y Q ) = f ty £ J*. (Vi (UjHo)) dh\ 
Jo 4 

where the right hand side is well-defined as Young intergral. 

Proof. J^Lq, DhU^L Q satisfy (at least jointly with U^l ) RDEs driven by X which 
allows, in essence, to use Lyons' limit theorem; this is discussed in detail in [2"r?ll27| . 
A little care is needed since the resulting vector fields are not bounded anymore. 
However we can rule out explosion and then localize the problem: the needed 
remark is that J£l q also satisfy a linear RDE of form 

dJ*_ = dM x (0* (y ) , t) ■ J t * o (y ) 

and explosion can be ruled out by direct iterative expansion and estimates of the 
Einstein sum as in [25^ . □ 



%n fact, G N (R d ) can realized as all points in the tensor algebra which arise from computing 
iterated integrals up to order N of smooth paths over a fixed time interval. The group product then 
corresponds to the concatenation of paths, the inverse corresponds to running a path backwards 
in time etc. 
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3. RDES DRIVEN BY GAUSSIAN SIGNALS 

We consider a continuous, centered Gaussian process X = [X , ,..,X d ) with 
independent components started at zero. This gives rise to an abstract Wiener 
space (W,W,m) where W = H C C* ([0,T] ,R d ). Note that H = ® d i=l H {{) and 
recall that element of Ti are of form h t = E (Y t £ (ft,)) where £ (ft.) is a Gaussian 
random variable. The ("reproducing kernel") Hilbert-structure on Ti is given by 
(ft,ft% :=E(e(ftK(ft')). 

Existence of a Gaussian geometric p-rough path above X is tantamount to the 
existence of certain Levy area integrals. The case of fractional Brownian motion 
is well understood and several construction have been carried out [7] [29] [11] [35] . 
In particular, one requires H > 1/4 for the existence of stochastic areas (which 
can be defined as L 2 (P)-limits as in Ito's theory). Resultingly, one has to deal 
with geometric p-rough paths for p < 4. (When p < 2 there is enough sample path 
regularity to use Young integration and we avoid speaking of rough paths.) 

Condition 1. Assume X lifts to a (random) geometric p-rough path X and 3q : 
1/p + l/q > 1 such that 

H^C q - var ([0,T},R d ) . 

The example to have in mind is Brownian motion for which the above condition 
is satisfied with p = 2 + e and q = 1. (We shall say more about other Gaussian 
examples in section [U) 

If X = R H denotes the geometric p-rough path, p e (l/-ff, [l/H] + 1), associated 
to fractional Brownian motion then it satisfies a Stroock- Varadhan support descrip- 
tion in rough path topology. This was first conjectued by Ledoux et al. [28] (who 
obtained it for the Brownian rough path) and carried out independently in [TT| [15] 
for H > 1/3. The case of H > 1/4 is more difficult and discussed in some detail 
in [4] . A proof in the generic context of Gaussian rough paths (covering fBM with 
H > 1/4 as special case) is given in the forthcoming paper .18] . The statement is 



(3.1) supp(P,X) = {S [p] (W)} 

where support and closure are relative to the homogenous p-variation topology for 
geometric p-rough paths. We recall that S[ p ], for [p] =2,3 given by 

#2 : h^l+ ( dh+ ( f dh(&dh 
Jo Jo Jo 

/>t />t />s />t />s />r 

S 3 : h^>l+ I dh + / dh<E>dh+ / dh<S)dh<E>dh 
Jo Jo Jo Jo Jo Jo 

lifts IR^-valued paths canonically to (M rf )-valued paths. In [18 it is seen that 
X exists provided the covariance has finite p-variation with p < 2 and it is also es- 
tablished that Ti. '—t C ,p_var which guarantees that S[ p ] (Tt) is well-defined via Young 
integration. Such support description will be important in checking condition [5] 
section 0] 

Definition 1. [5S] [35] [M] Given an abstract Wiener space (W 7 H,p), a r.v. F : 
W — ► R is H differentiable at u> g W iff exists DF (to) E H* such that 

VheH: i^-F{u} + eh)\ = (DF (cu) , h) n . 
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A vector-valued r.v. F = (F 1 , F e ) : W — ► M e is 7i differentiable iff each F l is TL 
differentiable. In this case, DF (lo) = [DF 1 (lo) , ...,DF e (lo)) is a linear bounded 
map from TL — > R e . One i/ien defines the Malliavin covariance matrix as the random 
matrix 

a(uj) := ((DF\DF%) i j=h e ER e * e . 
We call F weakly non- degenerate i/det (a) 7^ almost surely. 

Proposition 2. Assume condition {7]j. Then, for fixed t > 0, i/ie K e -TOZwe(f random 
variable 

^ ^ ^0 (yo) 

zs almost surely TL- differentiable. 

Proof. By assumption l/p+l/q > 1. We may assume that X (w) has been defined 
so that X (lo) is a geometric p-rough path for every uj £ W. We also know that 
X (lo + h) = T/jX almost surely and so 

P [X (lo + e n h) = X^X for any countable family (e n )] = 1 

We fix lo in the above set of full measure. For fixed t, define 

z itS = jH s (v t (u%_ )) . 

Noting that s <— > Z% s is in C p_var we have, with implicit summation over i, 



\D h U]l (y )\ = 



JlL. (Vi(uZ_ ))dh\ 



Zidti 



< c \Z 



\p—var 



< c 



Z(0)\)^\h\ p . var 
(\Z\ p - var + \Z(0)\)x\h\ H 



(DUH (y ),h) n 



(We used Young's inequality.) The linear map DU^_ (yo) : h 1— > DhU^L a (yo) £ 
is bounded and each component is an element of TL*, hence 

Noting that the derivative at e — exists, by definition, if the difference quotients 
converge as e J. and this holds iff convergence to the same limit takes place along 
any sequence e„ J, 0. It follows that, for almost every lo, 

d 



h ^ <! —U 1 

de 



rX(LJ+e^l) 



(yo)\ =(DU*L (y ),h) n 

> e=0 



and so the random variable U^_ (yo) is indeed a.s. 7i-differentiable. 



□ 



4. Conditions on Driving Process 

We now give a complete list of assumptions on the (d-dimensional) Gaussian 
driving signal (X t : t £ [0, T]). The first condition was already needed in the pre- 
vious section to show 7i-differentiability of RDE solutions driven by X; we repeat 
it for completeness and to give some additional examples. 
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Condition 2. Assume X lifts to a (random) geometric p-rough path X and 3q : 
l/p+ l/q > 1 such that 

(4.1) n^c g - var ([o,T],R d ) . 

In the Brownian motion case this holds, as already remarked earlier, withp = 2+ 
e and q — 1 . The same is true for the Brownian bridge and the Ornstein-Uhlenbeck 
examples discussed in the introduction; although case-by-case verifications are not 
difficult, there is general criterion on the covariance which implies (|4. 1 [) . see [17j 
Prop 16 applied with p = 1], which also covers fBM. Let us give a direct argument 
for case of fBM which covers any Hurst parameter H > 1/4. Writing TL H for the 
Cameron-Martin space of fBM, the variation embedding in [16] gives 

U H ^ for any g > (# + 1/ 2 ) ~ 1 . 

At the same time [3 [29l [TTJ [35] fBM lifts to a geometric p-rough path for p > 
1/H. By choosing p, q small enough 1/p + l/q can be made arbitrarily close to 
H + (H + 1/2) = 2H + 1/2 > 1 and so gl]) holds indeed for fBM with Hurst 
parameter H > 1/4. 

Condition 3. Fix T > 0. We assume non-degeneracy on [0,T] in the sense that 
for any smooth f — (/i, fd) ■ [0, T] — > M. d we have 




Again, fBM satisfies the following non-degeneracy condition simply because 
Cl ([0, T] , R d ) cH H ,ci. PH. A Brownian bridge which returns to zero at time T 
is ruled out, while a Brownian bridge which returns to zero after time T is allowed. 
Checking condition [3] for the Ornstein-Uhlenbeck example in the introduction is left 
as exercise for the reader. The following lemma taken from (4j contains a few ram- 
ifications concerning condition [3J since U = ef =1 W (i) there is no loss in generality 
in assuming d = 1. 

Lemma 1. Assume 

1/p + l/q > 1. 

(i) The requirement that f is smooth above can be relaxed to f € (jp~ var , 

(ii) The requirement that J fdh = OVTi £ Ji can be relaxed to the the quantifier "for 
all h in some orthonormal basis of TL" . 

(Hi) The non- degeneracy condition^ is equivalent to saying that for all smooth 
f =/= 0, the zero-mean Gaussian random variable J Q fdX (which exists as Young 
integral or via integration-by-parts) has positive definite variance. 

(iv) The non-degeneracy condition^ is equivalent to saying that for all times < 
ii < ... < t n < T the covariance matrix of (A tl , X tn ), that is, 

is (strictly) positive definite. 

(v) Non-degeneracy on [0,T] implies non- degeneracy on [0, t] for any t S (0, T]. 

Condition 4. "0-1 law": The germ a-algebra <~)t>ocr (X s : s £ [0,t]) contains only 
events of probability zero or one. 
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When X is Brownian motion, this is the well-known Blumcnthal zero-one law. 
More generally, it holds whenever X is an adapted functional of Brownian motion, 
including all examples (such as fBM) in which X has a Volterra presentation [5] 

X t = f K{t,s)dB s . 
Jo 

(Nothing is assumed on K other than having the above Wiener-Ito integral well- 
defined.) The 0-1 law also holds when X is the strong solution of an SDE driven 
by Brownian motion; this includes the Ornstein-Uhlenbeck - and Browian bridge 
examples. An example where the 0-1 law fails is given by the random-ray X : t 
tBx (ui) in which case the germ-event {u> : dX t (u>) /dt\t=o+ > 0} has probability 
1/2. (In fact, sample path differentiability at 0+ implies non-triviality of the germ 
CT-algebra see [10] and references therein) . We observe that the random ray example 
is (a) already ruled out by condition [3] and (b) should be ruled out anyway since it 
does not trigger to the bracket phenomenon needed for a Hormander statement. 

The next condition expresses some sort of scaled support statement at t = 0+ 
and is precisely what is needed in the last part (Step 4) in the proof of the Theorem 
[1] below. We give examples and easier-to-check conditions below. To state it, we 
recall 25, Thm 2.2.1] that a geometric p-rough path x lifts uniquely lifts uniquely 
and continuously (with respect to homogenous p-variation distances) to a path in 
the free step-iV nilpotent groupQ, say 

S N (x) e C£" var ([0, T] , G N (R d )) for N>\p\. 

We also recall that G N (R d ) carries a dilation operator S which generalizes scalar 
multiplication on K d . 

Condition 5. Assume there exists H G (0, 1) such that for all fixed N > \p\, 
writing X = S N (X), all g £ G N (R d ) and for all e > 0, 

liminfP (d ^ n flXi/„,sj < e) > 0. 

Proposition 3. Let B denote d-dimensional fractional Brownian motion with fixed 
Hurst parameter H S (1/4, 1) and consider the lift to a (random) geometric p-rough 
path, denoted by X = B, with p < 4. Then it satisfies condition^ 

Remark 1. Brownian motion is covered with H = 1/2. 

Proof. Write B = Sn (B). From section [3l and the references therein, the support 
of the law of B w.r.t. homogeneous p- variation distance is C °' p " var ([0, T] , (M (R d )) , 
that is, the closure of lifted smooth path started at with respect to homo- 
geneous p- variation distance [2S1[T7]. By continuity of Sn [HI Thm 2.2.1] fol- 
lowed by evaluation of the path at time 1 it follows that the support of the law 
of Bi is full, that is, equal to G N (R d ). On the other hand, fractional scaling 

(n H B t / n : t > 0) = (B t : t > 0) implies 5 n -EfBi/ n = Bi and so, thanks to full sup- 
port of Bi , 

liminfP (d (S nH B 1/n , g) < ej = P (d (Bi, 9 ) < ej > 0. 

□ 

^The in CQ -var indicates that Xq is started at the unit element in the group. 
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Although scaling was important in the previous proof, it is only used at times 
near 0+. One thus suspects that every other Gaussian signal X which scales sim- 
ilarly (on the level of N th iterated integrals!) also satisfies condition [5j To make 
this precise we need 

Theorem 2 ([17] ). Let (X, Y) = (X 1 , Y 1 , . . . , X d , Y d ) be a centered continuous 
Gaussian process on [0,1] such that [X 1 ,Y 1 ^ are independent for i = l,...,d. Let 
p G [1,2) and assume the covariance of (X,Y), as function on [0, l] 2 , is of finite 
p-variation (in 2D sens^). Then, for every p > 2p, X and Y can be lifted to 
geometric p-rough paths denoted X and Y. Moreover, there exist a constant C 
depending only onp,p, the covariance of (X,Y) so that for all q € [l,oo), 

\dp-var (X, Y)| i9 ^ < C^/q \Rx-Y\ oo;[0,1] 2 • 

(Note that Rx-y (s, t) is a diagonal matrix with entries depending on s, t.) 

Corollary 1. Let (X,B) satisfy the conditions of the previous theorem and assume 
that B is a (d- dimensional) fractional Brownian motion with fixed Hurst parameter 
H G (1/4, 1). Assume in addition that 

(4-2) ™ 2ff |^-sU[o,i/«]^0. 

Then condition^ holds. 

Proof. With focus on one diagonal entry and with mild abuse of notation (writing 
X,B instead of X\B\..) 

n 2H |-Rx-s|oo;[0,l/n] a 

= sup E[n H (X s/n - B s/n ) n H (X t/n - B t/n )\ 
s,te[o,i] 

which can be rewritten in terms of the rescaled process 

X (n) = n H X . /ni and simi- 
larly for B, as 



sup E 

s,te[o,i] 



= R 



;[o,i]' 



By assumption and the previous theorem, this entails that 
d p _ var (x (n) ,B (n) ) ^ in probability. 

By continuity of Sn, still writing X'™) = Sn (X^ n ') for fixed N, and similarly for 
BW, we have 

d (x[ n) ,B[ n) ^ < d p _ var;[0 ,i] (x (n) ,B (n) ) -> in probability. 



^Given a function / from [0, l] 2 into some normed space, its variation (in the 2D sense!) is an 
immediate generalization of the standard definition but based on "increments" of form 

A [a,b) x \c,d) = f(b,d) + f (a, c) - f(a, d) - f(b, c). 
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But then 



d[S n HX 1/n ,g) < e 



d(X?>,g)<e 



> Fid 



X.^,B^)+d(B^,g)<e) 



> 



d[B[ n \g)<s/2 



-P{d{X^\B ( ^) >e/2 



and so 



UmMF (d (s n nX 1/n , g^j < ej > liminfP (d (B^ ti) , < e/2 

and this is positive by the example in which we discussed the case of B resp. B. 
The proof is finished. □ 

Example 5 (Ornstein Uhlcnbcck). Given a Brownian motion B, we consider the 
Ornstein Uhlenbeck (short: OU) process given by the ltd integral 

X t = I e-^-^dBr 





If B is d- dimensional this yields, component-wise, the d-dimensional OU process. It 
is readily checked that (X,B) satisifes the assumptions of theorem^ (In fact, one 
sees p = 1 and we are dealing with geometric p-rough paths of Brownian regularity, 
i.e. p = 2 + e.) Condition H4-<ty then holds with H = 1/2 : take s,t £ [0, 1/n] and 
compute, with focus on one non-diagonal entry, 

Rx- B {s,t) = E[(X s -B s )(X t -B t )} 

(e-^ - l) (e~^ -l\dr = (n~ 3 ) . 

By corollary [7] we see condition [3| holds for the Ornstein Uhlenbeck examples. 

Example 6 (Brownian Bridge). Writing Xf := B t — t^Bt where B is standard 
Brownian motion this follows along the same lines, again by comparison with B t 
fort-> 0+. 

5. Taylor Expansions for Rough Differential Equations 

Given a smooth vector field W and smooth driving signal x (■) for the ODE 
dy = V (y) dx, it follows from (EH?)) that 

JgU (W {yf)) = W ( yo ) + f J *U ([Vi, W] (y:)) dx\, 

Jo 
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where Einstein's summation convention is used throughout. Iterated use of this 
leads to the Taylor expansion 



J^ t (W(yf)) = W\ V0 + [V U W] 



yo ~r " J lyo x o,t 
[V % ,[V 3 ,W]]\ yo ^ 

[V n ,...[V tN ,W]}\ yo . 



1;» 



where x , t denotes the signature of x (•) \L f , in R d © (M^)® 2 © • • • © (IR^)®^ © . . . 
Such an expansion makes immediate sense when x is replaced by a geometric p- 
rough pathj. Remainder estimates can be obtained via Euler-estimates |19j pro- 
vided Jq^ 4 (W (yf)) is a solution of some ODE of form dz = V (z) dx. This is 
accomplished by setting 

z := 0\ z\ z 3 ) := (y*, J %, J£_ t (W (yf ))) € IT © K exe © R e 

Noting that J§^ t (W (yf)) is given by z 2 ■ W (z 1 ) in terms of matrix multiplication 
we have 

dz 1 = Vi (z 1 ) dx i 

dz 2 = -z 2 ■ DVi (z 1 ) dx 1 

dz 3 = (dz 2 ) ■ W (z 1 ) + z 2 ■ d (W (z 1 )) 

= z 2 ■ (-DVi (z 1 ) ■ W (z 1 ) + DW (z 1 ) ■ Vi (z 1 )) dx i 

= z 2 -[V^WWzidx 1 

started from (yo, I, W (yo)) where / denotes the identity matrix in ]R exe and we see 
that V is given by 

Vi ( z 1 ) \ 
Vi (z\z 2 ,z 3 ) = I -z 2 -DV l (z 1 ) , i=l,...,d. 



z 



2 



[Vi,W] (zi) 



Lemma 2. Assume Vi, V4, W are smooth vector fields, bounded with all deriva- 
tives bounded. Then V — V^J is a collection of smooth (possibly unbounded) 
vector fields but explosion does not occur. More precisley, there exists a unique RDE 
solution to dz — V (z) dx on any compact time interval [0, T]. In fact, for some 
increasing function tp from M + into itself 

Moc;[0,t] < <P{M) When ll X llp-™r,-[0,t] < M - 

Proof. Smoothness of V is obvious and so the RDE dz = V (z) dx has a solution 
up to some possible explosion time. From the particular structure of V we now 
argue that explosion cannot occur in finite time: z 1 does not explode as it is a 
genuine RDE solution along Lip vector fields, z 2 does not explode as it satisfies a 
linear RDE (driven by some rough path M x as already remarked in the proof of 
Proposition [T|) . Clearly then, z 3 = z 2 ■ W (z 1 ) where W is a bounded vector fields 
cannot explode. More precisely, using the estimates for RDE solutions along Lip 
respectively linear vector fields in |19j respectively |25j it is clear that z remains in 



^By definition, such a p-rough path takes values in the step-[p] tensor algebra but recall that 
there is a unique lift to the step- TV group for any N > [p]. 



14 



THOMAS CASS AND PETER FRIZ 



a ball of radius only depending on M if ||x|| p var .j t , < M. (With some care one can 
show that logy (M) = 0(M P ) as M — > oo but this is irrelevant for the sequel.) □ 

Let us make the following definitions: given (to — l)-times differcntiable vector 



fields V — (Vi, V d ) on R e , g e 



3 fc=0 



and i/ 6 R e we write 



£(V)(y,g):=E E S k - lu -' lk V n ---V lk I(y). 

k=l 

G{l,...,d} 

(Here I denote the identitify function on R e and vector fields identified with first 
order differential operators.) In a similar spirit, given another sufficiently smooth 
vector field W we first write 

[V n ,V l2 ,---V lk ,W}:= [V n ,[V i2 ,--- [V ik , W] ■ ■ ■ ]] 

(which may be viewed as first order differential operator) and then 

(5.1) g k -[V,...,V,W]\ V0 := S k ' iu -' ik [Vn,V i2 ,---V ik ,W] 

ii_,...,i k 
e{i,...,d} 

with the convention that g° • T4 = Vk ■ 

Proposition 4. (Localized Euler Estimates) Given ||x|| p waT ..j t i < M and some 
integer to > p — 1 there exists C — C (M) = C = C (^M, m,p, V^j such that 

TTy (z , 0; x) o t - (zq, S m (x) o t J < C (M) x t^~ . 

Proof. If V were bounded with bounded derivatives this would be a consequence 
of [HI Thm 19]. On the other hand, z must remain in the ball -B (0, ip(M)) and 
we can replace V by (compactly supported) vector fields V such that V = V on 
B (0, tp (M )). After this localization we apply [19l Thm 19]. □ 

Lemma 3. Let f be a smooth function on R e lifted to a smooth function on R e © 

R exe 0R e by 

f(z\z 2 ,z 3 )=f(z a ). 
Viewing vector fields as first order differential operators, we have 

V il ---V iN \J=[V il ,--- ,V iN ,W]\ v J. 
As a consequence, if I denotes the identity function on R e , 



< 



z?-W\ yo -J2KflV,--- ,V,W] 

k=l 

z o,t — £(y) f z °' ( x )c 



Proof. Taylor expansion of the evolution ODE of z 3 (t) shows that Vi ± ■ ■ ■ Vi N \ Zo f 
[Vn,--- ,V iN ,W]\ yo f . 

Corollary 2. Fix a e T yo R e = R e with \a\ = l. Then 



□ 



"JSUt (W (yf)) -E« T (Kt ■ [V, ■ ■ ■ , V, W] 



k=0 



t=l/n 



-mH 



with n 
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Proof. We estimate this probability by 

m 

r JgU (W (y?)) - ]T « T (*o, t ■ [V, • • • , V, W] 



t=l/n 



e 

>-n" 



iff. 



l X llp-var;[0,l/n] — ^ 



l X llp-var;[0,l/n] — ^ 



then, using |a| = 1 and the previous lemma, 



< 



< 



z : l/n — £(y\ \^ Z 0, S m (x) 0j l/ n 

C(l)x 



>-n 



-mH. 



X llp-var;[0,l/?7.] — ^ 



>2 n 



o (1) using the localized Euler estimates. 



The probability of the (deterministic) event 

C(l) 









(a 




Vn) 









will be zero for n large enough provided m+1 > mi/ which is the case since p > 
l,i?<l. □ 

6. On Hormander's condition 

Let V — (Vi,...,Vd) denote a collection of smooth vector fields defined in a 
neighbourhood of yo G R e . Given a multi-index I — (ii, ...,ik) € {1, with 
length |/| = k, the vector field Vi is defined by iterated Lie brackets 

(6.1) ^ : =[y n ,i/ i2 ,...,y lfc ]^ [t/ n ,[t/ l2 ,...,[^_ 1 ,y i j...]. 

If W is another smooth vector field defined in a neighbourhood of yo € R e we writ^l 
^ ■ [V,...,V;W] := Yl a h '-' ik - 1 [V il ,V i2 ,...,V ik _ 1 ,W] 

e(M d )® (k - 1} leng'th fc 6{W} 

Recall that the step-r free nilpotent group with d generators, G"'(R d ), was realized 
as submanifold of the tensor algebra 

TM (R d )=®l =0 (R d f k . 
Definition 2. Given reN »e say £/ia£ condition (H) r holds at yo £ R e «/ 

(6.2) span{F/| yo : \I\ < r} - 7; j0 K e £* K e ; 
Similarly, we say that (HT) r holds at yo if the span of 

(6.3) J TTfc-x (g) • JV, ...,V, Vj] \ yo :k = 1, -, r; » = 1, rf, g e G^ 1 ^) 

^ length k 

equals T yo R e = K e - Hormander's condtion (H) is satisfied at yo iff (H) r holds for 
some reN. Similarly, we say that the Hormander-type condtion (HT) is satisfied 
at yo iff (HT) r holds for some reN. (When no confusion arises we omit reference 
to y Q .) 



®We introduced this notation already in the previous section, cf. ( 15.11 . 
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Proposition 5. For any fixed r £ N, the span of equals the span of W.3\) . 

Consequently, Hormander's condition (H) at yo is equivalent to the Hormander-type 
condtion (HT) at yo- 

Proof. Given a multi-index I = {i\, ik—i, ik) of length k < r and writing e±, ...ed 
for the canonical basis of R d 

g = g (*!,..., ifc-i) 

= exp (t^J <g> •• • ® exp (ife-ie^) 
€ G r_1 (R d ) C T r_1 (K d ) . 

(Recall that T r_1 (R d ) is a tensor algebra with multiplication ®, exp is defined 
by the usual series and the CBH formula shows that the so-defined g is indeed in 
G r ~ 1 (R d ) as claimed.) It follows that any 

TTfc-! (g)- [y,...,y,y<4 U 

length 

lies in the (HT) r -span i.e. the linear span of (|6 . 3[) . Now, the (HT) r -span is a closed 
linear subspace of T yo M. e = M. e and so it is clear that any element of form 

Kk-i(d a g)- [v,...,y,v ik ] \ yo 

length k 

where d a stands for any higher order partial derivative with respect to ti, 
i.e. 

9q =(-) withaeCNuW)*" 1 

is also in the (HT) r -span for any t±, . . . , tfc_i and, in particular, when evaluated at 
t\ = ■ ■ ■ = tk-i = 0. For the particular choice a = (1, . . . , 1) we have 

gk-l 



dt 1 ...dt k _ 1 



■g|ti=o....tfc_i=o - e-ii 



where h is an element of T' r_1 (R d ) with the only non-zero entry arising on the 
[k — l) th tensor level, i.e. 

7r fc _i (h) = e i± ® ••• ® e jfc _ 1 . 

Thus, 

TTfc-x (h) • jv; ... y, y it ]j TO = , . . ■ , v^,, , v ik ] \ yo 

length k 

is in our (HT) r - span. But this says precisely that, for any multi-index I of lenght 
k < r the bracket vector field evaluated at yo i.e. Vi\ yo is an element of our (HT) r - 
span. □ 

7. Proof of Main Result 

We are now in a position to give 

Proof (of Theorem^. We fix t € (0,T]. As usual it suffices to show a.s. invertibil- 
ity of 

a t = ((DY t \DY t j ) ) e K exe . 

V \ / Ti' i.j — 1 ,e 
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In terms of an ONB (h n ) of the Cameron Martin space we can write 



(7.1) a t = Y,( DY t' h ")H®( DY uK) H 



E 



f J* s (v k {¥„)) dh k niS ® / j* a w (n)) < ia 

o Jo 



(Summation over up-down indices is from here on tacitly assumed.) Invertibility of 
cr is equivalent to invertibility of the reduced covariance matrix 

Ct ~ E I J ^s (Vk (Y s )) <, s ® / J * s (y 8 )) <, s 
„ ./o 7o 

which has the advantage of being adapted, i.e. being a (X s : s g [0, immeasurable. 
We now assume that 

P (dot C t = 0) > 

and will see that this leads to a contradiction with Hormander's condition. 
Step 1: Let K s be the random subspace of T yo W = R e . spanned by 

{j£- r (V k (Y r ));re[0,s},k = l,...,d}. 

The subspace K + = r\ s> oK s is measurable with respect to the germ cr-algebra and 
by our "0-1 law" assumption, deterministic with probability one. A random time 
is defined by 

9 = inf {s e (0,t] : dimK s > dimi^ +} A t, 
and we note that 6 > a.s. For any vector v e M e we have 

r* 2 

v T C t v = T / v T J^ s (Vk(Y s ))dh k n , s . 

Assuming v T Ctv — implies 

Vn: / v T J^ s (Vk(Y s ))dh k niS = 
Jo 

and hence, by our non-degeneracy condition on the Gaussian process 

v T J^ s (Vk (Y s )) = 

for any s € [0,i] and any k — l,...,d which implies that v is orthogonal to K t . 
Therefore, K + ^ M. e , otherwise K s = R e for every s > so that v must be zero, 
which implies C t is invertible a.s. in contradiction with our hypothesis. 
Step 2: We saw that K + is a deterministic and linear subspace of R e with strict 
inclusion K + C R e i n particular, there exists a deterministic vector z e M e \ {0} 
which is orthogonal to K Q + . We will show that z is orthogonal to to all vector 
fields and (suitable) brackets evaluated at yo, thereby contradicting the fact that 
our vector fields satisfy Hormander's condition. By definition of 9, K + = K t for 
< t < and so for every k = 1, ...d, 

(7.2) z T J^ t (V k (Y t )) = for t < 9. 

Observe that, by evaluation at t = 0, this implies z _L span{Vi, Vj} \ yo . 

Step 3: We call an element g € ®™ =0 (M. d )® k group-like iff for any N e N, 



(tto (g) , . . . , n N (g)) G G N (R d ) C ©£Lo W 
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We now keep k fixed and make induction hypothesis I (m — 1) : 

Vg group-like, j < m - 1 : z T vr, (g) [V, ■ ■ ■ , V; V k ] \ yo = 0. 

To this end, take the shortest path 7™ : [0, 1/n] — ► M d such that S m (7™) equals 
7r i,...,m (g), the projection of g to the free stcp-m nilpotent group with d generators, 
denoted G m (R d ) . Then 

l7"li. var;[ o,i/„] = lki,...,™(g)|| Gm(Rd) < 00 

and the scaled path 

h n (t) = n~ H j n (t) , He (0, 1) 
has length (over the interval [0, 1/n]) proportional to n~ H which tends to as 
n — > 00. Our plan is to show that 

(7.3) Ve > : limmfP ( | z T j£_ 1/n (v k (y>?] n ) ) | < e/n mH ) > 

which, since the event involved is deterministic, really says that 

n mH z T Jti /n (v k (y$ n ))\<e 

holds true for all n > no (e) large enough. Then, sending n — > 00, a Taylor expan- 
sion and I (m — 1) shows that the l.h.s. converges to 



z n 7r m (S m (h n )) ■ [V, •• ■ ,V; V k ]\ 



= 7Tm(g) 



< e 



and since e > is arbitrary we showed I (m) which completes the induction step. 
Step 4: The only thing left to show is (|7.3|) . that is, positivity of liminf of 



> 



v{\^j£i/n{Vk(lfi n ))\<e/Tr B ) 

P (jz T J^. (V k (y.)) - z T 4i. (v k (/")) I =i/n < e/n^ 



-P(6 < 1/n) 
and since 6 > a.s. it is enough to show that 

z T J6^.(v k (y-))-z T 4Z. (v k (y h ' 



lim inf P 

n — >oo 



•=l/n 



< e/n 



mH 



> 0. 



Using I(m— 1) + stochastic Taylor expansion (more precisely, corollary [5]) this is 
equivalent to show positivity of lim inf of 

p (|z T x™. [v,---,V; v k ] - z T 4i. (y k (v h "))\ =i/n < s/n mH ^ . 

Rewriting things, we need to show positivity of liminf of 



P(\n m »z+ [V,--- ,V;V k ]-K™ 1/r , 
or, equivalently, 



z T n m 



H Jot 1/n (V k (y$ n ))\ 



-^z T [V,-,V;V k ]K m (g) 



P (\z T [V, ■ ■ ■ , V; V k ]\ yg (n mH X™ 1/n - n m (g)) | < e) 
But this is implied by condition [5] and so the proof is finished. 



□ 
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